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e xt t n dt for 

Large Values of n ] 

Walter Gautschi 

An asymptotic expansion is given which is well suited for numerical computation when 
n is large and x arbitrary positive. 

1. Let 

/» 00 

E n (x)=\ e- xt t~ n dt; z>0; n=l, 2, 3, . . . . (1) 

By means of four integrations by parts, G. Blanch 2 has found the approximation 

w x _^! |" u n | n(n—2x) n(Sx 2 — 8nx+n 2 ) ~\ 

* nKX) ~x+n\_ ^(z+n) 2 (x+ri) 4 ' "*" (x+^) 6 J [ } 

She also gives an integral representation for the error. Formula (2) has proved very efficient for 
computing E n (x) for large values of n in the whole range x>0. In what follows, the eomplete 
expansion is given, as well as error estimates. 

Denote by h k (u) the polynomial (of degree Jc-1 if &>0), defined recursively by 

h k+l (u) = (\-2ku)h k (u) + u(l+u)h k {u) (Jc=0, 1,2, ... ), h (u)=l. (3) 

Let 

and let <x k , fi k be lower and upper bounds, respectively, for H k (u) in the interval u>0: 

* k <H k (u)<p k (u>0). (5) 

Then it will be proved that 

^ )= £[§M£^~' + **H' (6) 

a*n-* < E k (x,n) < fc (i+-^JL_)>-\ (7) 



1 This paper was prepared under a National Bureau of Standards contract with The American University. 

/» 00 

2 G. Blanch, An asymptotic expansion for E n (x)= I (e- xu ju n )du, NBS Applied Math. Series 37, 61 (1954). 
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For reference, the first eight polynomials h k (u) and corresponding values of a k , fi k are listed : 3 

h (u) = h 1 (u)=l 
h 2 (u) = l — 2u 
hz(u)=\-8u+<5u 2 
K(u)=l-22u+58u 2 —2Au d 
h 5 (u)=l-52u+328u 2 -U4:U z +120u 4 
h 6 (u)=l-U4:U+U52u 2 -U00uz + l 3708u i -720u 5 
h 7 (u) = l-24:Ou+5610u 2 -32120u"+58UOu 4 -^9S4u 5 +504:Ou & 
ai = 0, a 2 =-0.07, a 3 =-0.18, a 4 = — 0.36, a b =— 0.60, a 6 =-0.94, a 7 =-\A 

/3l = & = & = ft=ft=ft=l, ^7=1.8. 

2. Consider, more generally, the integral 

/= T e fW dt, (8) 

U a 

where f(t) is a real function defined on the finite or infinite interval (a, b). It is assumed that 
J(t) has derivatives of any order in (a, b) and that f (t)^0. Following van der Corput and 
Franklin, 4 we define the sequence g k (t) by 

9o(t)= T(f)' gk+l{t) =ff) {k=0 ' 1 ' 2 ' ■ ■ - } - (9) 

Setting 

h=C g k (t)f(t)e^dt, 

clearly Iq=I, and integration by parts yields 

I*=MtW w ]l- ( b g'ic(t)e f «\n=v k -I k+1 , (10) 

J a 

where 

v k =g k (b)e m -g k (a)e^. (11) 

Hence, 

I=v -v l +v 2 -v 3 + . . . +(-iy-%- l +(-iyi k , (12) 

I*= F gl-i{t)S m dt. (13) 

J a 



3 The author is indebted to Mrs. L. K. Cherwinski and Mrs. B. H. Walter for the calculation of the otk and Pk. 

* J. G. van der Corput and Joel Franklin, Approximation of integrals by integration by parts, Nederl. Akad. Wetensch. Proc. Ser. [A] 54 
213-219 (1951). 
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In case of an infinite interval {a, b) it has to be assumed that the values (11) exist for all k. 
Equation (10) then shows that the existence of the integral I k implies the existence of I k+X . 

3. The integral in (8) is equal to E n {x) if 

f(t) = — (xt+nlxit), a=l, 6—00. 

A short computation shows that with this definition of f(t), the sequence g k (t) in (9) is 
equal to 

,_. (-l) k+ H h k (u) _, (-l) k+ H „,. _ k xt 

xt+n (l+u) 21 * xt-f-n n 

where the h k (u) are the polynomials defined in (3) and H k (u) the rational functions (4). For the 
quantities v k in (11), one obtains 



v k = 
Furthermore, 






^i_ l «) = (-l)*H»(it)n-* > «-$■ 

Hence, from (12) and (13), 

^ )= £[§K^^ (H) 

By (5) 

* k E n (x) <rH k (^\e- xt r n dt<p k E n (x) f 

and using the well-known inequality, 5 

<e x E n (x)< * t (*><», 



X + 'tt £ + 71— 1 

one gets 






*+ 






From this and (14) the formulas (6), (7) follow immediately. 

It may be observed that the result (6), (7) holds also for nonintegral values of n with n>l. 



5 E. Hopf, Mathematical problems of radiative equilibrium, Cambridge Tracts in Mathematics and Mathematical Physics, No. 31, p. 26 
(Cambridge University Press, 1934). 



Washington, October 14, 1958. 
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